Quantum dissociation of an edge of a Luttinger liquid 
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In a Luttinger liquid phase of one-dimensional molecular matter the strength of zero-point mo- 
tion can be characterized by dimensionless De Boer's number quantifying the interplay of quantum 
fluctuations and two-body interactions. Selecting the latter in the Morse form we show that dissoci- 
ation of the Luttinger liquid is a process initiated at the system edge. The latter becomes unstable 
against quantum fluctuations at a value of De Boer's number which is smaller than that of the bulk 
instability which parallels the classical phenomenon of surface melting. 
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A classical three-dimensional soHd melts through a 
first-order transition at a temperature when the free en- 
ergies of the solid and liquid phases coincide. At suffi- 
ciently low temperatures quantum effects dominate and 
a quantum solid can melt due to zero-point motion 0. 

The most curious feature of classical melting is the dif- 
ficulty in overheating the solid while supercooling the liq- 
uid is easy. The latter is expected for the first-order tran- 
sition while the former is explained by the phenomenon 
of surface melting: often, as the bulk transition is ap- 
proached, the melting begins at the free surface of a solid. 
The surface melting is well-documented experimentally, 
and phenomenologically it can be viewed as a wetting of 
the solid by its own melt Q. 

A well-understood example of surface melting of a 
quantum solid is that of the edge melting of the two- 
dimensional Wigner crystal in a strong magnetic field 0| 

The goal of this note is to point out that strictly one- 
dimensional matter with a free edge can also exhibit an 
analog of surface melting. Fundamentally this happens 
because the edge represents a zero-dimensional system 
subject to stronger quantum fiuctuations than the one- 
dimensional bulk. Due to broken translational symmetry, 
zero-point motion modifies the cohesive properties of the 
edge differently from those of the bulk. 

Experimentally one-dimensional matter can be real- 
ized in carbon nanotube bundles 0|. The latter can 
play a role of one-dimensional hosts for foreign particles 
that can find themselves bound in the interstitial chan- 
nels or inside the tubes Qj. Additionally one-dimensional 
atomic chains can be constructed on selected templates 
with the help of scanning tunneling microscopy, or via 
self-assembly of the deposited material of the chain . 

Consider a many-body system of identical particles of 
mass m with pairwise interaction V(h) corresponding to 
the molecular matter 01 : at large interparticle separation 
h the interaction is dominated by weak rapidly decaying 
van der Waals attraction, while at short distances there is 
a strong overlap repulsion . As a result, the pair potential 
V(h) has an asymmetric minimum at some intermediate 
h. Assume that the pair potential is of the form 



where e is the energy scale of the potential, I is the po- 
tential range, Qq is a dimensionless parameter, and U{y) 
is a function common to a family of substances. With 
this choice the quantum theorem of the corresponding 
states fl| holds stating that every property measured in 
appropriate dimensionless units is only determined by the 
function U{y), particle statistics and De Boer's number 
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measuring the intensity of zero-point motion. 

The possibility of several bulk phases in the system 
translates into the corresponding number of the branches 
of the energy as a function of Aq; the lowest of them 
singles out the ground state of the system. When two 
energy curves cross, the ground state changes via a first- 
order phase transition. For sufficiently large Aq and zero 
pressure the ground state must correspond to individual 
particles infinitely far apart from each other. This is a 
monoatomic gas which will be chosen as the zero refer- 
ence point for the energy. 

In what follows we select the pair interaction potential 
in the Morse form lojl: 



V{h) = e(e"2(V;~Qo) _ 
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V{h)^eU{h/l-Qo), 
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where e is the depth of the potential well and Qo is 
the location of the minimum of ijSl measured in units of 
the potential range I. Similar to the applications of the 
Lennard-Jones potential to laboratory molecular systems 
01, the only reason behind this choice is the possibility 
of analytic progress. Morse parameters for a series of 
molecular substances and corresponding De Boer's num- 
bers were computed in Ref.|l(l|. Hereafter the energy 
and length scales will be measured in units of e and I, re- 
spectively. As appropriate for molecular substances, we 
restrict ourselves to nearest-neighbor interactions. 

In the classical limit, Aq = 0, the ground state of the 
system is a crystal; its quantum counterpart for suffi- 
ciently small Ao is a Luttinger liquid whose properties 
have been computed in Ref.0| as follows: 

The length of any bulk bond ft. as a function of imag- 
inary time T is viewed as a quantum-mechanical degree 
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of freedom subject to the external potential V{h). This 
bond joins together two half-infinite segments represent- 
ing the rest of the system, the "bath". After the bath 
is approximated by a harmonic liquid, the latter can be 
integrated out away from the anharmonic bond leading 
to a problem of the Caldeira-Leggett type The lat- 
ter has been analyzed by a combination of variational 
and renormalization- group techniques, and it has been 
demonstrated that the approximation is a controlled way 
of dealing with the interplay of zero-point motion and 
anharmonicity of the two-body interaction Similar 
consideration applied to the edge bond of a half-infinite 
Luttinger liquid leads to the EucHdian action of the form 
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where p and c are the mass density and sound velocity, 
respectively, and h{uj) is the Fourier transform of the 
bond-length field; the frequency cutoff is given by the 
Debye frequency ujd- 

The pc|w| form of the kinetic energy term of the action 
Q can be understood by noticing that if the bond length 
oscillates with frequency w, then during one oscillation 
period 2n/\uj\ this disturbance propagates in the bulk a 
distance of order c/\uj\. Thus the usual kinetic energy 
density, proportional to puj^ should be multiplied by the 
size of the region c/|w| affected by the motion. 

The calculation of the properties of the edge of a Lut- 
tinger Hquid proceeds through the appHcation to the ac- 
tion Q of Feynman's variational principle 0| which 
states that for any trial action 5*0 with associated ground- 
state energy Eq, the system's true ground-state energy is 
bounded above by Eq + (T/h) < S — Sq >o where the 
zero-temperature limit T = is taken at the end and 
<>o denotes an expectation value computed with So- 

Similar to the bulk problem the trial action is se- 
lected in the Gaussian form 



pc 
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duj 
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where dimensionless variational parameters Q and 7 have 
a meaning of the bond length and its stiffness, respec- 
tively. Then the root-mean-square (rms) fiuctuation of 
the bond length can be computed as 



where 



A = 



TTpCp 



(6) 
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quantifies the strength of zero-point motion in the Lut- 
tinger liquid. The binding energy of the edge particle 



Eedge is approximated by Eq + {T/h) < S — Sq >q, i. e. 
by its upper bound: 

Eedaeh,Q) = (7rAg/A)ln(l + 7)-2e««"'?(l+7-')'^ 
+ e2(Q«-Q)(l +7-1)8^ (8) 

Minimizing E^dge with respect to Q we arrive at the ex- 
pression for the quantum expansion of the edge bond 

Qedge-Qo = 6Aln(l+7-i) (9) 

Substituting this back into ijHJ, E^dge can be written as 

Eedgeh) = M^/A) ln(l + 7) ^ (1 + 7^')^'^ (10) 

Minimizing En. lfT(1|l with respect to 7, and substituting 
the outcome back into ifTTHl we find 



(11) 



and 



E.dgc = (^A2/4A2) (4A ln(l + 7) - 7) (12) 



respectively. The results ©, ®, CHj and lEJ should 
be compared with their bulk counterparts ,1(1] : 
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?6uifc-(3o = 3Aln(l + 7r/2), 



Ao = A(l + 7r/2)- 



(1 + 7r/2)-2^ (ttA ln(l + 2/7r) - 1) (16) 
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Substituting En. lfT5|l back in Ens. ltTT|l and lfT2|l brings 
them into a form convenient for analysis 



7=(4A)(l + 7r/2)2^(l + 7-^) 



, -1N-4A 



(17) 



Eedge - (7r/4)(l + 7:/2)-'^ (4A ln(l + 7) - 7) (18) 

The properties of the edge as a function of the quantum 
parameter A {JJ can be computed by finding a solution 
7(A) to Ea. (ll7|l minimizing the energy Ijl8|l and substi- 
tuting the outcome in the expressions for the rms fiuc- 
tuation ^ and quantum expansion (EJ; the dependence 
on De Boer's number JJl follows from En. lll5|l . 

In the classical Hmit, A ^ 0, the only solution to itTzjl is 
7 — A/tt with the energy III8II Eedge = — 1 as expected. As 
the degree of zero-point motion intensifies (A increases) , 
the bond stiffness 7 decreases and the energy Eedge in- 
creases. For finite A En. ltTTjl may have more than one 
solution. One of them is always 7 = corresponding to 
the delocalized edge particle. For large A this solution 
must correspond to the lowest (zero) energy II18|I . 
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For 7 <C 1 the right-hand-side of Ea. (jl7ll behaves as 
7^'^ while for 7 ^ oo it approaches a 7-independent limit, 
thus implying that ifTzIl cannot have more than three so- 
lutions and that A = 1/4 plays a special role. 

For < A < 1/4 Eq.(Q3l has two solutions and the 
larger of them (whose A = limit is 7 = 4/7r) corre- 
sponds to the lowest energy (tlSll . For A = 1/4 the explicit 
solution to is 7 = (4/7r)(l + - 1 ~ 1.0415 

As A increases beyond 1/4, Eq.^T7| acquires a third 
root whose A 1/4 H- limit is 7 = {Tr/4{1 + 
7r/2)i/2)i/(4A-i) _^ Q However this solution leads to 
a larger energy itTsll than even the delocalized solution 
7 = 0. The lowest energy (bound) state continues to be 
described by the largest solution to Ijl7|l . 

As A continues to increase, the finite solutions to ifTzjl 
approach each other and at some A they coalesce. This 
is a critical phenomenon corresponding to the limit of 
stability of the bound edge. At that point the slopes of 
the right- and left-hand-sides of En. ltTzjl coincide which 
leads to the limiting values 7 ~ 0.4920 and A ~ 0.3730 
satisfying the relationship 7 = 4A — 1. At larger values of 
A Ea. lll7|l has only one solution 7 = corresponding to 
an unbound edge. The transition between the bound and 
unbound states actually happens before the limit of sta- 
bility is reached, namely when the energy ifTsjl vanishes. 
Numerical analysis shows that it happens at A ~ 0.3412. 
This is close to the limit of stability thus implying that 
the edge delocalization is a weak first-order transition. 

The results of the analysis are summarized in Fig. 1 
where we show the bond stiffness 7 and the edge bind- 
ing energy Eedge as functions of the quantum parameter 
A. The metastability develops in the 1/4 < A < 0.3730 
range: for A < 0.3412 the bound edge has lower energy 
while for A > 0.3412 the ground state corresponds to a 
delocalized edge particle. 

These conclusions should be contrasted with the prop- 
erties of the bulk Luttinger liquid. Its range of exis- 
tence is given by < A < 1.0591 (or equivalently 
< Ao < 0.3896) which is the condition that a solution 
A(Ao) to Ea. lll5ll can be found for given De Boer's num- 
ber Ao ©. Therefore in the 0.3730 < A < 1.0591 range 
the bulk Luttinger liquid is stable against the disorder- 
ing effect of quantum fluctuations while the edge is not. 
This is due to the stronger softening effect that zero-point 
motion has on the free edge as compared to the bulk of 
the system. The direct evidence of this is presented in 
Fig. 2 where we show the quantum expansion and rms 
fluctuation of the bulk and edge bonds as functions of 
the quantum parameter A 10 within their corresponding 
ranges of existence. The quantum expansion in the bulk 
H14II is a linear function of A while the edge bond expands 
faster than linearly because the bond stiffness 7 entering 
the argument of the logarithm in (jsj is a decreasing func- 
tion of A as shown in Fig. 1. Since the 7(A) dependence 
is not very strong one can say that the edge expansion is 
roughly twice the bulk value as suggested by the ratio of 
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Figure 1: Dimensionless stiffness of the edge bond 7 and 
corresponding binding energy of the edge particle E^dge of 
a half-infinite Luttinger liquid as functions of the quantum 
parameter A iQ. The region of metastability is confined to 
the 1/4 < A < 0.3730 range. 
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Figure 2: (Color online) Quantum expansion Q ~Qo and rms 
fiuctuation of the bulk and edge bonds as functions of the 
quantum parameter A Q. The rms fluctuation is shown both 
as the vertical extent of shaded regions centered around the 
quantum expansion curves, and explicitly in the inset. 



pre-logarithmic factors in Eqs.ljHJ and H14II . This can be 
understood by noticing that any bulk bond joins two half- 
inflnite Luttinger liquids thus implying that its dynamics 
is twice as inertial as that of the edge. In this sense zero- 
point motion at the edge is about twice as strong as that 
in the bulk. The same argument explains why the edge 
rms fluctuation is roughly square-root of two larger than 
its bulk counterpart (compare Eqs.JHJ and Ill3p l. 

In describing the dynamics of the edge bond the rest of 
the system was approximated by a harmonic liquid with 
the bulk properties which means that the bond adjacent 
to the edge has the length and rms fluctuation identical 
to those in the bulk. This is an artifact and in reality, as 
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Figure 3: The dependences of the energy per particle for var- 
ious bulk phases of the system on De Boer's number Ao Q 
together with edge binding energy. The arrow pointing down 
is the dimer dissociation threshold. 

one goes inside the bulk, the bond lengths and their rms 
fluctuations decrease approaching the bulk values asymp- 
totically. This deflciency would be acceptable provided 
the calculated length of the edge bond and its rms fluctu- 
ation are not very different from their bulk counterparts. 
Since for molecular matter with pair interaction potential 
of the Morse form the classical bond length satisfies the 
condition Qo ^ 5 lO], inspection of Fig. 2 shows that 
even at the limit of its stability the length of the edge 
bond and its rms fiuctuation do not exceed their bulk 
counterparts by more than an acceptable 25%. 

Moreover, the relative fluctuation < p >o^edge I Qedge 
is always signiflcantly smaller than unity which implies 
that our conclusions are weakly sensitive to the statis- 
tics of the underlying particles and that the deficiencies 
of the Morse potential in mimicking the true pair in- 
teraction at largest and shortest distances are ignorable. 
The latter allows us to argue that the edge dissociation 
pre-emting the bulk instability is a general property of 
one-dimensional molecular matter. 

In order to gain an insight into the consequences of this 
effect in Fig. 3 we plot the ground-state energy per par- 
ticle of the bulk Luttinger liquid (given by Eas. ljl5|l and 
Ijlfill l and the binding energy of the edge particle (de- 
termined through Ens. ltTfill . ifTTji . and itTsjl as functions 
of De Boer's number Ao Additionally we show the 
ground-state energy per particle for an infinitely diluted 
gas of Morse dimers, i?dimer(Ao) = -(l/2)(l-7rAo/V2)^ 
The bold parts of the curves describe the ground 
states of the bulk matter: as De Boer's number increases, 
at Ao — 0.3365 the Luttinger liquid evaporates via a 
discontinuous transition into a gas of dimers followed 
by a continuous dissociation transition at Ao — yf^ji^ 
into a monoatomic gas For a system with a free 

edge the binding energy of the edge particle i?edge can 
become smaller than its dimer counterpart Edimer- for 



Ao ^ 0.1981 the whole Luttinger liquid comes unraveled, 
two particles at a time despite the fact that the bulk con- 
densed state is energetically favorable. Since our bulk 
and edge binding energies are variational upper bounds, 
in actuality the dimer gas may not come into play; its 
role then will be played by the monoatomic gas. 

H the escape of the edge particles to infinity is impos- 
sible due to a distant obstacle, this will generate a vapor 
pressure and the bulk Luttinger liquid may coexist with 
a gas of particles. As Aq increases toward the point of the 
bulk transition, dissociation proceeds inside the bulk in a 
manner similar to that in surface melting Q . We hasten 
to mention the speculative character of the statements of 
this paragraph which we plan to clarify in the future. 

The examples of one-dimensional matter with dissoci- 
ated edge and stable bulk, 0.1981 < Ao < 0.3365, include 
H2 and D2 in free space, and more cases can be found in 
the presence of a medium • 
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